Acceleration of heavy and light particles in turbulence: comparison between 
experiments and direct numerical simulations 



R. Volk/'H E. Calzavarini, 2 G. Verhille, 1 D. Lohse, 2 N. Mordant, 3 J.-F. Pinton, 1 and F. Toschi 4 

1 Laboratoire de Physique, de I'Ecole normale superieure de Lyon, 

CNRS UMR5672, 46 Allee d'ltalie, 69007 Lyon, France 
^Faculty of Science, J. M. Burgers Centre for Fluid Dynamics, 
and Impact-Institute, University of Twente, 7500 AE Enschede, The Netherlands 
8 Laboratoire de Physique Statistique de I'Ecole normale superieure de Paris, 
CNRS UMR8550, 24 rue Lhomond, 75005 Pans, France 
Jf Istituto per le Applicazioni del Calcolo CNR, Viale del Policlimco 137, 00161 Roma, Italy 
and IN FN, Sezione di Ferrara, Via G. Saragat 1, 1-44 1 00 Ferrara, Italy. 

We compare experimental data and numerical simulations for the dynamics of inertial parti- 
cles with finite density in turbulence. In the experiment, bubbles and solid particles are optically 
tracked in a turbulent flow of water using an Extended Laser Doppler Velocimetry technique. The 
probability density functions (PDF) of particle accelerations and their auto-correlation in time are 
computed. Numerical results are obtained from a direct numerical simulation in which a suspension 
of passive pointwise particles is tracked, with the same finite density and the same response time 
as in the experiment. We observe a good agreement for both the variance of acceleration and the 
autocorrelation timescale of the dynamics; small discrepancies on the shape of the acceleration PDF 
are observed. We discuss the effects induced by the finite size of the particles, not taken into account 
in the present numerical simulations. 

PACS numbers: 47.27. Jv,47.27.Gs,02.50.-r 



I. INTRODUCTION 

Understanding the transport of inertial particles with 
finite density, such as sediments, neutrally buoyant par- 
ticles or bubbles in turbulent flows of water is of practical 
interest for both industrial engineering or environmental 
problems. In a turbulent flow, the mismatch in density 
between the particles and the fluid causes light particles 
to be trapped in high vortical regions while heavy parti- 
cles are ejected form vortex cores and concentrate in high 
strain regions As particles with different buoyancy 
tend to concentrate in different regions of the flow, they 
are expected to exhibit different dynamical behaviors. In 
recent years, significant progress has been made in the 
limit of infinitely heavy, pointwise particles 0, Q, and 
numerical simulations have received experimental sup- 
port [J, |5j. In case of infinitely light particles (bubbles): 
the result of the numerical simulations onparticlc dis- 
tributions and on fluid velocity spectra 0, 0, El agree in 
various aspects with experimental findings [9L TloL flU [HJ 
although direct comparison between experiments and nu- 
merical simulations for the acceleration pdf and correla- 
tion of the particles as not been investigated in the past. 

Indeed, in spite of the growing resolution of Direct 
Numerical Simulations (DNS) of the Navier-Stokes equa- 
tions at high Reynolds numbers, it remains a challenge to 
resolve the motion of realistic inertial particles: some de- 
gree of modelization is necessary. The equation of motion 
of finite-size, finite-density particles moving in a turbu- 
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lent flow is not precisely known, and a comparison with 
experimental data can help refining the models and ex- 
tending their range of validity. 

Several experimental techniques have been developed 
for measuring the velocity of particles along their tra- 
jectories. The optical tracking method developed in the 
Cornell group has revealed that fluid particles experience 
extremely intense accelerations (l3j , while individual par- 
ticles have been tracked for time durations of the order 
of the flow integral time scale using an acoustic tech- 
nique [T4[ . Because of the very fast decrease of the acous- 
tic scattering cross-section with the scatterer's size, this 
method is limited to particles with diameter of the order 
of the wavelength, i. e. inertial range sizes [To] , |23| . The 
principle of the acoustic technique is completely analo- 
gous to laser Doppler velocimetry (LDV), provided that 
expanded light beams are used (an arrangement we call 
E-LDV hereafter). The advantage of E-LDV, compared 
to acoustics, is that the much smaller wavelength of light 
allows a better resolution in space and also the use of 
smaller tracer particles. Theprinciple of the measure- 
ment technique is reported in [16[ , where its performance 
has been compared and validated against the silicon-strip 
tracking (T3 , Il7j of neutrally buoyant Lagrangian trac- 
ers. We focus here on the dynamics of inertial particles 
i.e. particles whose density differ from that of the fluid. 
Here we report the first comparison between experimen- 
tal measurements of acceleration of particles having a 
relative density in the range 1(T 3 (air bubbles) to 1.4 
(PMMA) in the same highly turbulent flow, and numer- 
ical results obtained by tracking pointwise particles with 
finite density in a direct numerical simulation of isotropic 
homogeneous turbulence [l||, Qj| . 



2 




e ~ 25 W/kg 
R\ ~ 850 



FIG. 1: Experimental setup. (Top left): schematics of the von 
Karman flow in water - side view. (Top right): principle of 
the Laser-Doppler Velocimetry using wide beams (ELDV) - 
top view of the experiment. PM: location of the photmultiplcr 
which detect scattering light modulation as a particle crosses 
the interference pattern created at the intersection of the laser 
beams. 



Numerical simulations are performed by means of stan- 
dard pseudo-spectral methods, where particular care has 
been used in keeping a good resolution at the dissipative 
scales. The numerical code for integrating the evolution 
of the Eulerian field and the Lagrangian tracing of par- 
ticles is the same as described in 0, H, H3 ■ A thorough 
validation of the numerical approach, included the La- 
grangian evolution of the tracers has recently been per- 
formed against experimental measurements [261 ] . The nu- 
merical integration of tracers has, with respect to exper- 
iments, the clear advantage of a uniform, well controlled 
geometry and very large statistics; on the opposite, the 
resolution can be limited to small Reynolds numbers. 
For what concerns the treatment of realistic particles, 
i.e. particles with a density mismatch and a "finite" size, 
the best modelization to use is not clear and one of the 
main goal of this manuscript is indeed to compare state- 
of-the-art Lagrangian data against numerical results from 
a current modelization. 



II. EXPERIMENTAL SETUP AND RESULTS 

The Laser Doppler technique is based on the same 
principle as the ultrasound Doppler method which has 
good tracking performance of individual Lagrangian trac- 
ers [3, [HI. In order to access dissipative scales, and in 
particular for acceleration measurements, one adapts the 
technique from ultrasound to laser light: the gain is of a 
factor 1000 in wavelength so that one expects to detect 
micron-sized particles. For a Lagrangian measurement, 
one has to be able to follow the particle motion to get 
information about its dynamics in time. For this, wide 
Laser beams arc needed to illuminate the particle on a 
significant fraction of its path. The optical setup is an 



extension of the well known laser Doppler velocimetry 
technique; Fig. [1] A Laser beam is split into two beams; 
each is then expanded by a telescope so that their diam- 
eter is about 5mm. Then the two beams intersect in the 
flow where they create an array of interference fringes. As 
a particle crosses the fringes, the scattered light is modu- 
lated at a frequency directly proportional to the compo- 
nent of the velocity perpendicular to the fringes. It yields 
a measurement of one component of the particle velocity. 
In practice, we use a CW YAG laser of wave length 532 
nm with 1.2 W maximum output power. In order to get 
the sign of the velocity we use acousto-optic modulators 
(AOM) to shift the frequency of the beams so that the 
fringes are actually travelling at a constant speed. The 
angle of the two beams is tuned to impose a 60 microns 
interfringe so that the frequency shift between the beams 
(100 kHz) corresponds to 6 m/s. As the beams are not 
focused, the interfringe remains constant across the mea- 
surement volume whose size is about 5 x 5 x 10 mm 3 . It 
is imaged on a photomultiplier whose output is recorded 
using a National Instrument PXI-NI5621 digitizer at rate 
1 MHz. 

The flow is of the Von Karman kind as in several pre- 
vious experiments using acoustics [14j or optical tech- 
niques [l3j . Water fills a cylindrical container of internal 
diameter 15 cm, length 20 cm. It is driven by two disks 
of diameter 10 cm, fitted with blades in order to increase 
steering. The rotation rate is fixed at values up to 10 Hz. 
For the measurements reported here, the Taylor based 
Reynolds number reaches up to 850 at a maximum dis- 
sipation rate e equal to 25 W/kg. We study three types 
of particles: neutrally buoyant polystyrene particles with 
size 31 microns and density 1.06, PMMA particles with 
size 43 microns and density 1.4 and air bubbles with a 
size of about 150 microns. The mean size of the bubbles, 
measured optically by imaging the measurement volume 
on a CCD, is imposed by the balance between the inter- 
facial surface tension a and the turbulent fluctuations of 
pressure. This fragmentation process is known to lead to 
a well defined and stationary size distribution [20| with 

a typical diameter D cx (cr/p/) 3 ^ 5 e~ 2 / 5 , pj being the 
density of the fluid. 

The signal processing step is crucial as both time and 
frequency - (i.e. velocity) - resolutions rely on its per- 
formance. Frequency demodulation is achieved using the 
same algorithm as in the acoustic Doppler technique. It 
is a approximated maximum likelihood method coupled 
with a Kalman filter [2l| : a parametric estimator assumes 
that the signal is made of a modulated complex exponen- 
tial and Gaussian noise. The amplitude of the recorded 
signal and the modulation frequencies are assumed to be 
slowly evolving compared to the duration of the time win- 
dow used to estimate the instantaneous frequency. Here 
the time window is about 30 ^s long and sets the time 
resolution of the algorithm. Outputs of the algorithm 
are the instantaneous frequency, the amplitude of the 
modulation and a confidence estimate which is used to 
eliminate unreliable detections. Afterwards, the accclcr- 
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Heavy 
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1 
1 

0.75 
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0.31 
4.1 
1.03 
1.64 


2.85 ±0.07 
2.94 ±0.07 
2.63 ±0.12 
25.9 ±0.46 


1 

1.03 
0.92 
9.08 



TABLE I: (top) Parameters of the particles in the von 
Karman flow at R x = 850 (r/ = (u 3 /e) 1/4 = 17 /im and 
t v — \Jv /e — 0.26 10~ 3 s). p p and p/ are the densities of the 
particles and fluid, and r p = a 2 /(30u) is the stokes response 
time of the particles. The Taylor-based turbulent Reynolds 
number is computed as R\ = ^ 'I5uf ms / oj measuring the 
one-component root-mean-square velocity, it rma , with the E- 
LDV system and e by monitoring the power consumption of 
the motors. The non dimensional constant ao is derived from 
the Heisenberg-Yaglom relationship. The last column com- 
pares the value for the inertial particle to the one obtained 
for the Lagrangian tracer (which is denoted by the subscript 
T). (bottom) Same as above: parameters of the particles 
tracked in the DNS of homogeneous isotropic turbulence at 
R\ = 180. Out of the numerically analyzed 64 parameter 
combinations (0,St), we have picked those which were close 
to the experimental values for (/3, St). 



ation of the particle is computed by differentiation of the 
velocity output. Note that measurements are performed 
only when a particle moves within the (limited) measure- 
ment volume so that after processing, the data consists 
in a collection of sequences with variable lengths. For all 
the measurements, the acceleration variance is computed 
using the same procedure as in [l7j : it is obtained for 
several width of the smoothing kernel used in the differ- 
entiation of the velocity signal and then interpolated to 
zero filter width. 

For small neutrally buoyant particles, i.e. Lagrangian 
tracers, our data is in excellent agreement with the high 
speed im aging m easurements performed by the Cornell 
group 

liallSllll- When the rms value of the acceleration 
is normalized by the Heisenberg-Yaglom scaling: (a 2 ) = 
(e beeing the energy dissipation rate per unit 
mass and v = 1.3 • 1CP 6 m 2 s _1 the kinematic viscosity 
of the fluid), both experiments yields the same values 
for the non dimensional constant ao (ao = 6.4 ± 1 at 
R x = 850 for the E-LDV compared to 6.2 ± 0.4 for the 
Cornell data at R x = 690). 

We have applied our technique to compare the dy- 
namics of Lagrangian tracers to the one of heavier or 
lighter particles (sec table U for numbers). We first com- 
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FIG. 2: Probability distribution function of accelerations, 
normalized to the variance of the data sets, (top) Data from 
experiment at R\ = 850. (middle) DNS of homogeneous 
isotropic turbulence at R\ = 180. (bottom) Comparison of 
experimental measurements and DNS results. 



pute the velocity root mean square value u rms for the 
three cases: the values are {1.1, 1.2, 1.0} ± O.lrn.s" 1 at 
R\ = 850 for the tracers, heavy (PMMA spheres), and 
light particles (bubbles). Within error bars, the large 
scale dynamics seems to be unaffected by changes in the 
particle inertia. The acceleration distribution and auto- 
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correlation in the three cases are shown respectively in 
Fig. [U (top) and Fig. [3] (top). The acceleration PDFs 
are quite similar for moderate acceleration values (below 
about 10a rms ), as also observed in low Reynolds number 
numerical simulations [221 ]. However, the probability of 
very large accelerations seems to be reduced in the case 
of inertial particles as compared to Lagrangian tracers. 
The normalized acceleration variance ao varies very sig- 
nificantly: it is reduced to 4.3 ± 1 for heavier particles 
while it is increased to 26 ± 5 for bubbles. The correla- 
tion functions also show significant changes with inertia: 
the characteristic time of decay is longer for heavy par- 
ticles and shorter for bubbles compared to tracers. We 
measure T COII /T n = {0.5, 0.9, 0.25} respectively for trac- 
ers, heavy, and light particles, with the correlation time 
defined as the half- width at mid amplitude of the correla- 
tion function. We thus observe important changes in the 
dynamics, even if the distribution of acceleration weakly 
changes with inertia. 

Note that in our setup the Kolmogorov length is about 
i] = 17 /im at R\ = 850, so that the bubble size is about 
10 r\ and therefore may not be considered as a point par- 
ticle. Thus, one may wonder if the bubble dynamics is 
not altered by spatial filtering as recently demonstrated 
for particles with diameters in the inertia range [23|. To 
check, we have compared the dynamics of large neutrally 
buoyant particles with diameter 250 fim to the one of 
Lagrangian tracers. The results is shown in Fig. [2] to- 
gether with the other particles: the effect of the parti- 
cle size on the PDF is found to be weak as the curve 
nicely superimposes with the ones for inertial particles. 
However, the size effect is clear when comparing either 
the coefficient ao (reduced to 2.2), or the autocorrelation 
functions. One observes that the correlation time of the 
large particles is twice the one for the tracers. We con- 
clude that the bubbles size may have a leading effect on 
the acceleration variance, and that the value of ao re- 
ported here probably underestimates the one that would 
be measured for smaller bubbles (with diameters closer 
to the Kolmogorov scale). 



III. COMPARISON WITH NUMERICAL 
SIMULATIONS 

We compare the experimental data with the results 
from a direct numerical simulation [l8|, [l9[ where a pas- 
sive suspension of pointwise particles with finite density 
are tracked in a homogeneous isotropic turbulent flow. 
The dynamics of the particles is computed in the most 
simplified form of the equation of motion, i.e. assum- 
ing that the particles are spherical, non-deformable and 
smaller than the Kolmogorov length scale of the flow [24| • 
When one retains only the Stokes drag force and the 
added mass effect, the equation of motion then reads: 

dv p Du 1 



where v p = x(t) is the particle velocity, u(x(t), t) the ve- 
locity of the fluid at the location of the particle described 
by the Navier-Stokes equation, while [3 — 3pj/(pf + 2p p ) 
accounts for the added mass effect and and t p = a 2 / (3f3v) 
is the Stokes response time for a particle of radius a. 
When made dimensionless by the Kolmogorov dissipa- 
tive scales (r r; , it^) eq. ([1]) reads 

dv c „D\i 1 , , . . 

with the particle acceleration a now expressed in the 
Heisenbcrg-Yaglom units. Thus, at a given Reynolds 
number, the particles dynamics only depends on the val- 
ues of the two dimensionless parameters (3 and St = T p /r rj 
This is generally different from the case of infinite inertia 
of the particles ((3 = 0), which has been formerly ad- 
dressed in several numerical and theoretical studies 0], 
and for which instead only the Stokes number St matters. 
It is also different from the pure bubble case {(3 = 3) for 
which theparticle indeed has no inertia but only added 
mass [E S Hi • We performed numerical simulation at 
Re\ = 180 (grid resolution 512 3 ), in which many parti- 
cles, characterized by different pairs (/?, St) (specifically 
64 different sets of O(10 5 ) particles) were numerically in- 
tegrated by means of cqn. ((T|). Particles do not have 
feedback on the flow field. 

In order to compare the numerical results with the ex- 
perimental data, three types of particles (tracers, heavy, 
and bubbles) with different inertia and Stokes number 
have been studied. The values for both (3 and St have 
been chosen close to the ones of the particles used for the 
E-LDV (see table HI . The evolution of the normalized ac- 
celeration variance shows the same trend in experiments 
and numerics: ao is reduced from the tracer value 2.85 to 
2.63 for heavier particles and increased to 26 for bubbles 
(table IJ). This seems to be a robust trend in the DNS. 
To emphasize this, in figure 0] we show the behavior of 
y/a~o, i.e. the root-mean-square value of the particle ac- 
celeration normalized by the Heisenberg-Yaglom scaling, 
in a wide range of the (/3, St) parameter space from a 
less turbulent DNS (Re\ = 75) which has a very large 
number of (/?, St) pairs. Results from the Re\ = 180, 
not showed here, arc qualitatively similar. Note again 
that no significant Reynolds number dependence of the 
probability distribution was found in ref. [l6j |. 

The acceleration distribution behavior and its compar- 
ison with the experiment is reported in Fig. [2j In the 
numerics we observe that the probability of very large 
accelerations is reduced for the heavier particles as com- 
pared with tracers, while it is increased for the bubbles. 
This feature, seems not to be present in the experimental 
results. Furthermore, we notice that for the three types of 
particles, the acceleration PDFs, rescaled by the rms ac- 
celeration, is close to the experiments. Experimental ones 
have always longer tails, reflecting the more intermittent 
nature of the turbulent flow, which has a larger Reynolds 
number (Re\.EXP = 850 vs. Rc\,dns = 180 ). We also 
observe a qualitative agreement for the changes in the ac- 
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FIG. 3: Autocorrelation coefficients of the accelerations: 
(top) Data from experiments at R\ = 850. (bottom) Data 
from DNS of homogeneous isotropic turbulence at R\ = 180 
For the (j3, Si) values we refer to table [I] 



celeration autocorrelation functions when changing iner- 
tia, figurc[3j One measures t cott /t v = {0.95, 1.35, 0.25} 
respectively for tracers, heavy and light particles. Just 
as observed for the experiments, the dynamics is faster 
for the bubbles while heavier particles decorrelate slower 
than fluid tracers. The Re\ difference is more pro- 
nounced here than in the pdf's (see ref. [HI) and prevent 
from a more detailed comparison. 



IV. DISCUSSION 

While solving a simplified version for the equation of 
motion, the numerics reproduces qualitatively the effect 
of the particles inertia on their dynamics. In particu- 
lar, the dependence of the acceleration autocorrelation 
on the particle inertia is nicely reproduced, see figure [3] 
However, also some discrepancies become visible, though 
not yet completely conclusive, as a better resolution and 
statistics of both the experiments and the numerics would 
be important for firmer conclusions. Nevertheless, in this 
section we shall have a closer look at the differences and 




FIG. 4: Behaviour of the normalized root mean square accel- 
eration y/a^ = ((a 2 )e" 3/2 ^ 1/2 ) 1/2 as a function of both St and 
f3 for a Re\ — 75 DNS. Iso-contour for v / ao7r (red) and the 
line y/ao,T ■ P + const, (green) are also reported. Note that ao 
does not depend on St for neutral (/3 = 1) particles. While 
it is always reduced/enhanced for heavy/light particles. For 
large particles (St ~ 4.1) we find ^/oq ~ pja^f. 



propose some explanations. 

First of all, there is only qualitative agreement on the 
ratio ao,ff/ao,T- It is larger for the experiment than for 
the numerics. Moreover, the tails of the numerical PDF 
of the bubble acceleration seem to be enhanced as com- 
pared to those for tracer acceleration. Vice versa, the 
tails of the numerical PDF of the particle acceleration 
seem to be reduced as compared to those for tracer ac- 
celeration. 

What is the origin of the difference between the ex- 
periments and the numerics? First of all the Taylor- 
Reynolds numbers are different, but ref. [l6| suggests an 
at most weak dependence of the acceleration PDFs on 
the Reynolds number; a finding that is supported by a 
comparison of our numerical simulations at Re\ = 185 
and Re\ — 75. 

Next, in the numerical simulations we disregarded the 
lift and the gravitational force. While this presumably 
has little effect on heavy particles and tracer, it does 
modify the dynamics of the bubbles. In refs. @, HI we 
had shown by comparison of numerical simulations for 
point bubbles with and without lift, that without lift 
the bubble accumulation inside the vortices is more pro- 
nounced, i.e., bubbles without lift are more exposed to 
the small scale fluctuations, which clearly will contribute 
to the pronounced tails of the numerically found acceler- 
ation PDF, see figurc[31 bottom. 

Next, also the two-way coupling of the particles (i.e., 
the back-reaction of the particles on the flow due to their 
buoyancy difference) has been neglected in the simula- 
tions of this paper. As e.g. shown in refs. 0, H| for bub- 
bles and in ref. [27| for particles, it has an effect on the 
turbulent energy spectrum and thus also on the acceler- 
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ation statistics. However, as in the present experiments 
the particle and bubble concentrations are very low, the 
two-way coupling effect on the spectra should be hardly 
detectable. 

The final difference between numerics and experiments 
we will discuss here - and presumably the most relevant 
one - is the finite size of the particles in the experiments 
as compared to the numerics which is based on effective 
forces on a point particle. Although the heavy particles 
are not large as compared to r\, this clearly holds for the 
bubbles and the 250/^to diameter neutral particles. In- 
deed, figure [3] shows how the finite size of these particles 
smears out the acceleration autocorrelation, as compared 
to the tracer case. Also the ratio ao,7v/ao,T for large neu- 
tral particles is only 0.34, which demonstrates that the 
size of large particles has a large effect on their acceler- 
ation variance. This type of spatial filtering, which also 
lowers the PDF of large neutral particles in the experi- 
ment, is not related to a temporal filtering of the particle 
based on its response time. This is clearly visible in Fig. 
[^middle) where one can see that two neutral particles 
(/3 = 1) with different response times (different St or r p ) 
have the same acceleration PDF, with same ao, and same 
autocorrelation function. Thus this size effect, which is 
not taken into account in the point-particle based numer- 
ical simulations, presumably is responsible for both the 
relatively small value of ao,s/oo,T measured for bubbles, 
and the change in the shape of the PDF. 

To conclude, we have reported acceleration measure- 



ments of inertial particles using extended Laser Dopplcr 
velocimetry and have compared the experimental data to 
DNS simulations of the motion of pointwise particles with 
finite density. We have observed a qualitative agreement 
between experiments and numerics in the shape of the 
PDF and of the autocorrelation function. We have given 
arguments for the small discrepancies. An experimen- 
tal study of the motion of bubbles with smaller sizes is 
needed for a better comparison with the numerical simu- 
lations. Also numerical simulations keeping into account 
the finite size of particles would presumably improve the 
agreement between experiments and numerical data and 
detailed comparisions as the one presented in this pa- 
per help to reveal the limitations of point-particle mod- 
els. Obviously, going beyond point-particles is extremely 
challenging in numerical simulations. A first step in this 
direction has e.g. been taken by Prosperetti and cowork- 
ers with their Physalis method [28l which presently is 
extended towards turbulent flows [29( . 
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